Methods of calculation of a friction coefficient: 
Application to the nanotubes 
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In this work we develop theoretical and numerical methods of calculation of a dynamic friction co- 
efhcient. The theoretical method is based on an adiabatic approximation which allows us to express 
the dynamic friction coefBcient in terms of the time integral of the autocorrelation function of the 
force between both sliding objects. The motion of the objects and the autocorrelation function can 
be numerically calculated by molecular-dynamics simulations. We have successfully applied these 
methods to the evaluation of the dynamic friction coefBcient of the relative motion of two concentric 
carbon nanotubes. The dynamic friction coefficient is shown to increase with the temperature. 

PACS numbers: 61.48.+C, 62.20.Qp 



Thanks to recent developments in nanotechnology, the 
hope is high to build mechanical devices on the scale 
of the nanometer. For this purpose, it is important to 
determine the mechanical properties and, especially, the 
friction forces in such nanodevices. In this work, we are 
going to focus on systems of carbon nanotubes. Our par- 
ticular interest is for a system that Cumings and Zettl 
observed experimentally with a TEM. They fixed an edge 
of a multiwalled nanotube to a surface and opened the 
other edge, they then extracted inner layers from the core 
for several nanometers and released them. They observed 
a full retraction of the inner layers and furthermore, they 
could conclude that the multiwalled nanotubes are self 
cleaning since amorphous carbon due to the opening of 
the edge are inside the tube and also do not have any 
wear, structural change or fatigue after several extrac- 
tion and retraction processes. This experiment suggest us 
that the multiwalled nanotubes can be promising systems 
for future nanometric mechanical parts such as springs, 
gears or even motors. 

A short time after the work of Cumings and Zettl, 
Zheng and Jiang Q estimated the frequency of the os- 
cillations in this system to be of the order of GHz. This 
result is also very interesting since in the macroscopic 
world moving parts with such frequencies does not ex- 
ist at the present time. But as in the macroscopic world 
moving parts have friction forces which hinder the motion 
and dissipate energy. We hence have to know the impor- 
tance of these forces before conceiving such devices. Our 
work will thus focus on, firstly methods of calculation of 
the friction and secondly an application of these theories 
to the multiwalled nanotubes. 

The plan of the paper will be as follows, we are first 
going to introduce the theoretical framework for the de- 
scription of the mechanics of nanotubes. We will then 
solve the classical equations of motion and calculate the 
dynamic friction coefficient by molecular-dynamics sim- 
ulation and by the autocorrelation-function method de- 
veloped by Jarzynski, Berry and Robbins 0, 

We have depicted in Fig. n]two shding nanotubes. R 
is the distance between the centers of mass of each nan- 



otube. The Hamiltonian of the system of two nanotubes 




FIG. 1: Sliding nanotubes. 



is given by 
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where Ti and T2 are respectively the total kinetic energies 



of the inner and outer tubes, V^g and V^g are the total 
Tersoff-Brenner potentials of both tubes. We use the set 
of parameters which optimizes the geometric structure 
0. The last term is the interaction potential between 
the tubes, which is taken as a 6-12 Lennard- Jones poten- 
tial calculated by Lu and Yang 4] and widely used for 
molecular-dynamics simulations of nanotubes |3j S IM 01 ■ 
The calculation of the friction coefficient will be done 
by using an adiabatic approximation since there are two 
different time scales in the system. The fast time scale 
is the one of the vibrations of the atoms which is of the 
order of the fs and the slow time scale is the one of the 
motion of the mass centers of the tubes which is of the 
order of the ps in our case. Thanks to this difference in 
time scales, we can use a method developed by Jarzynski, 
Berry and Robbins 8, 9] which is based on the fact that 
the following quantity is an adiabatic invariant [Tol [Tl| , 



.(2) 



n{E,et)^ J dz e[E-H{z,et)] 
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where e is a small parameter, E the energy of the fast 
system and H is an ergodic adiabatic Hamiltonian. This 
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means that H slowly evolves with time and that, at 
any instantaneous time, H produces trajectories which 
explore their energy shells ergodically and chaotically, 
z = (q, p) represents the fast phase-space variables. The 
smallncss of the parameter e guarantees that the mo- 
tion of vibration is fast compared to the motion of the 
mass centers. The other important assumption is that 
the force on the slow system can be expressed as an aver- 
age over the fast degrees of freedom since the time scale 
of the fast system is much smaller. We hence write the 
force as 

F ^ - J dz pVH (3) 

where p is the probability density of the fast system and 
V = d/dR is the derivative with respect to the slow 
variable. In our specific case R is the distance between 
the mass centers of the tubes. Using the fact that there 
are two different time scales, p is expanded as 

p(z, t) = po (z, R) + epi (z, R, t) + • • • (4) 

where the first order term po depends only on the slow 
time which is regarded as the parameter R. The evolu- 
tion of p is given by the Liouville equation. Injecting the 
expansion |0J in the Liouville equation gives a hierarchy 
of equations which can be solved using a microcanonical 
distribution for po- Having calculated p to the first order 
in e one can find the following reaction force, 

F = -\/E - r • R (5) 

where the first order term is the Born-Oppenheimer force 
with the energy of the fast system E and the second term 
the friction force where the velocity R of order e appears. 
The dynamic friction coefficients T are given by 

Ty- ^Q^dE (^dEn dr a, (r)^ (6) 

where Se is the partial derivative with respect to the en- 
ergy of the degrees of freedom of vibration E and Vl{E, R) 
is defined by Eq. ((2Jl. The autocorrelation function is 
given by the following expression 

C,,{t) - {d,HrdjH)E.,n (7) 

where H = H{z, R) — £'(R) represents the fluctuations of 
the Hamiltonian. We can simplify the expression of the 
friction in Eq. by using the Boltzmann equation S = 
fee IuOe^ and the thermodynamic relation dE/dS = T. 
Moreover, the adiabatic parameter R is one-dimensional 
in our system since the relative motion of the tubes is 
mainly in the longitudinal direction. We hence get 

r={P + dE) / dr C{t) ^13 dT C{t) (8) 
Jo Jo 

with P = l/k-QT. In Eq. we have used the result 
that the second term with Se is a negligible correction of 



the order of the inverse of the system size. The autocor- 
relation function is given by 

C(t) = {F^3{t)FU^))e,r ~ {Fui)Im (9) 

where i^Lj is the total Lennard- Jones force in the trans- 
verse direction, 
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We observe that Eq. ((SJ satisfies the fluctuation- 
dissipation theorem relating the fluctuating force acting 
on a system in a bath to the friction kernel in the gener- 
alized Langevin equation: 

A:bT C(i) = (/(t)/(0)) (11) 

with / — Flj — (i^Lj)- Eq. © is hence obtained in the 
Markovian limit where the friction kernel Hll|l decays on 
a short time scale such that C(t) = CoH^)- 

We are now going to show the results of the molecular- 
dynamics simulation. We focus our numerical studies on 
a specific nanotube namely (5,0)@(15,0) which is com- 
mensurate in the classification of Ref. 7] . The inner tube 
has iVi = 60 atoms and a length of = 10. 9A. The outer 
tube has N2 = 240 atoms and a length of I2 = 15. 3 A. 
The relative mass is hence p = 576 amu. The molecular- 
dynamics simulation is straightforward. We first extract 
the internal nanotube such as the initial distance between 
both mass centers is Rq — 6.0A. The tubes centers of 
masses have vanishing initial velocity and angular mo- 
mentum. We then start the simulation and solve the 
Newtonian dynamics with a fourth-order Runge-Kutta 
integration scheme in the microcanonical ensemble. The 
initial temperature of the system is 300K, which is intro- 
duced by a Maxwell-Boltzmann distribution of the initial 
velocities and a short early time relaxation. We use a 
time step of O.lfs. 

First, we can predict the energy to be dissipated by 
the sliding motion and thus the increase of the tem- 
perature of the system. Calculating the difference of 
energies between the minimum of the total Lennard- 
Jones potential Vl,j(-R) and the initial value VLj(i?o) 
gives us AEs = 0.623eV. Using AE = SNksT where 
N = Ni + N2 we hence should have an increase in tem- 
perature of AT = 8.047^. 

We observe in Fig. |21a double exponential decay in the 
motion. The origin of the double exponential decay holds 
in the fact that, in the initial configuration, a part of the 
inner tube is outside the outer tube. In this case, the 
friction is higher because the inner tube interacts with the 
edge of the outer tube and bounces during the reentry. It 
is only after a while that the inner tube oscillates inside 
the outer tube. We see that the oscillations are damped 
as expected. We can use a simple model to describe the 
dynamics when the inner tube remains inside, in which 
case the potential energy of the relative motion between 
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FIG. 2: (a) Evolution of the distance between mass centers 
represented by the successive maxima at each oscillation, (b) 
Evolution of the energy of the relative motion also repre- 
sented by the successive maxima at each oscillation. In this 
small system, the oscillations have a period of about lOps. 
The crossover between the exponentials happens around 120fs 
when the bouncing of the edges ceases. The inner tube is en- 
tirely inside the outer one after about 400fs. 

both tubes can be approximated by a harmonic potential. 
We hence have a Newton equation of the form, 

^iR = -TR-kR (12) 

which can be solved to show that the successive maxima 
in position R and energy E ~ i/ii?^ + ^kR? decay as 

|i?|max cx e-"/2M and i;,,ax«e-"/^ (13) 

We hence expect to have a ratio of 2 between the damp- 
ing rates of the maxima in energy and position. By per- 
forming fits of the curves in Fig. [21 we get respectively 
0.00072/ps and 0.0015/ps, in good agreement with the 
expectation. Multiplying by the relative mass /i we fi- 
nally get the friction coefficient V = 0.84 ± 0.02 amu/ps. 

The next step is to calculate the friction coefficient 
with the autocorrelation-function method. According to 
Eqs. (|S|)- (|10|I . the friction coefficient is given by the time 
integral of the autocorrelation function of the force acting 
on the slow system with the distance between the centers 
of mass fixed as a parameter. The calculation of the 
autocorrelation function is also performed by molecular- 
dynamics simulation except that the distance between 



FIG. 3: (a) Autocorrelation function of the intertube force, 
(b) Time integral of the autocorrelation function (a) numeri- 
cally converging to the friction coefficient. The small fluctua- 
tions seen beyond about 80fs in the autocorrelation-function 
(a) and its integral (b) give a negligible contribution and are 
reduced if the statistics used in the averaging are increased. 

the centers of mass is kept constant by compensating at 
each time step the motion of each atom by the motion of 
the mass center of their respective nanotube. We suppose 
that our system is ergodic and mixing (meaning that its 
autocorrelation functions vanish asymptotically at long 
times). Since we have nearly one thousand degrees of 
freedom this should be a good approximation. We can 
thus use the ergodic theorem and replace the ensemble 
average in Eq. Q by a time average. 

Figure O depicts the results of the numerical calcu- 
lation of the autocorrelation-function and its time inte- 
gral giving asymptotically the friction coefficient. With 
this autocorrelation-function method, we find the value 
r = 0.82 amu/ps for the friction coefficient, in excellent 
agreement with the value obtained by direct molecular- 
dynamics simulation. The great advantage of the auto- 
correlation function method is that the calculation is now 
about 50 times faster. Furthermore the autocorrelation 
function method is much more precise since it does not 
involve fitting to get the friction coefficient. 

Furthermore, the autocorrelation-function method al- 
lows us to obtain a theoretical formula for the friction 
coefficient. We first notice that the autocorrelation func- 
tion decays on the time scale of the vibrations of the 
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atoms in the carbon nanotubes. This time scale is of 
the order of the inverse of the Debye frequency wd so 
that the correlation time should be given approximately 
by tc ~ 2t: /luy). Since the Debye temperature of the 
nanotube is 8d ~ lOOOK 12], we hence expect the au- 
tocorrelation function to decay to zero over a correlation 
time approximately equal to 50fs, which is indeed con- 
firmed by Fig. 13 Moreover, the autocorrelation function 
at time zero is equal to the variance AF^ of the force. 
We can therefore approximate the time integral of the 
autocorrelation function and obtained the following ap- 
proximate formula for the friction coefficient: 



1 1 27r , , 



k^T 2 WD 



(14) 



Thanks to this formula, the evaluation of the friction 
coefficient reduces to the evaluation of the variance of 
the fluctuating force. On this ground, we now investigate 
the dependence of the friction coefficient on temperature. 
We have depicted in Fig. 0] the result of the numerical 
calculation of the standard deviation of the force. We see 



FIG. 4: Standard deviation of the force versus temperature. 

that the standard deviation of the force increases linearly 
with the temperature, hence the friction coefficient also 
increases linearly with temperature: T T . The fact 
that friction increases with temperature in our system 
can be understand by the fact that the dissipation of 



energy of the slow motion is due to the coupling with the 
vibrating degrees of freedom. Increasing the temperature 
results into larger vibrational motion and thus into higher 
friction. Furthermore, the calculation shows that edge 
effects give to the small system we here considered, a 
larger friction than to a system of double length. Then 
the friction increases in systems from double to longer 
lengths, the coefficient remaining of the same order of 
magnitude as in the small system here studied. 

In this work, we have developed theoretical and nu- 
merical methods to calculate the dynamic friction co- 
efficient in the one-dimensional sliding motion of two 
concentric carbon nanotubes. We focused our numeri- 
cal studies on a commensurate double-walled nanotube. 
We obtained by a molecular-dynamics simulation the en- 
ergy damping and thus the increase in temperature of the 
system. The simulation shows the importance of edge ef- 
fects at early times: a significantly larger damping is due 
to bouncing effects when the inner nanotube reenters. 
This effect disappears at long times when the inner tube 
moves inside the outer tube. The long time behavior is 
in good agreement with the one-dimensional model we 
proposed. We have then calculated the friction with the 
autocorrelation-function method. We observed that the 
autocorrelation function decreases to zero with a time of 
the order of 50fs which is a value in agreement with the 
experimental value of the Debye temperature. The inte- 
gral of the autocorrelation function gives the same result 
as the direct molecular-dynamics simulation with much 
less calculation time, actually a factor 50, and with bet- 
ter precision since we avoid the use of fitting methods. 
We can thus conclude that the autocorrelation-function 
method is very efficient for the calculation of friction co- 
efficients. Moreover, we have obtained a formula which 
relates the dynamic friction coefficient to the Debye fre- 
quency and the standard deviation of the fluctuations of 
the forces between the nanotubes and determined that 
the dynamical friction coefficient increases with the tem- 
perature. 
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